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Abstract. We construct initial data for a particular class of Brill wave metrics 
using Regge calculus, and compare the results to a corresponding continuum 
solution, finding excellent agreement. We then search for trapped surfaces in 
both sets of initial data, and provide an independent verification of the existence 
of an apparent horizon once a critical gravitational wave amplitude is passed. Our 
estimate of this critical value, using both the Regge and continuum solutions, 
supports other recent findings. 
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1. Introduction 



Regge calculus |l| is used to re-examine a particular class of Brill wave spacetimes 
investigated by Miyama || , and more recently, Holz et al M and Alcubierre et al 




In particular, we study the formation of apparent horizons within the simplicial initial 
data, and confirm recent findings regarding the appearance of such horizons (^] . 

Brill wave initial data was first constructed using Regge calculus by Dubai H, 
although we show elsewhere |jj that his procedure is unable to capture the full 
structure of complex axisymmctric initial data sets. By using an improved lattice 
(which is described in full in Ref. 0), we are able to accurately reproduce the 
continuum solution, and find apparent horizons in excellent agreement with recent 
two and three-dimensional calculations |qj. 

We proceed as follows. In the next section, we briefly survey the approach taken 
by Brill [|| in the construction of gravitational wave initial data. In section [| we 
describe the axisymmetric lattice used to construct Brill waves using Regge calculus. 
In section |] we describe the method used to find apparent horizons at a moment of 
time symmetry, and in section ^ we investigate the appearance of apparent horizons 
in the Brill wave space as the wave amplitude is increased. 
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a 


Mr 


M c 


M (Holz et al ) 


M (Alcubierre et al ) 


1 


3.38 x 10" 2 


3.38 x 10" 2 


3.40 x 10" 2 


3.38 x 10- z 


2 


1.26 x 1CT 1 


1.26 x KT 1 


1.26 x lO" 1 


1.27 x 10" 1 


5 


6.98 x 1CT 1 


6.98 x 10- 1 


6.96 x 10" 1 


7.00 x 10- 1 


10 


2.91 


2.91 


2.91 


2.91 


12 


4.67 


4.67 


4.67 


4.68 



Table 1. Mass estimates for the Regge (Mr) and continuum solutions (Mq). 
The mass is calculated from the decay of ip in the asvmptotic region. We also 
show the results of previous calculations by Holz et al g and Alcubierre et al Q , 
which are in excellent agreement. All our results were calculated on a 600 x 600 
grid. 



2. The continuum model 

We examine the class of axisymmetric vacuum spacetimes containing gravitational 
waves which admit a moment of time symmetry originally studied by Brill Q . Brill 
relaxed the conformal flatness criteria and introduced a spatial metric of the form 

ds 2 = ip i {e 2q (dp 2 + dz 2 )+p 2 dcb 2 }, (1) 

where the arbitrary function q(p, z) can be considered the distribution of gravitational 
wave amplitude In order to obtain an axisymmetric solution which is reflection 
symmetric about z — 0, we require q. p = along z = and q = q_ z = on p = 0, 
together with the condition that q falls off sufficiently quickly (at least r -2 ) in the 
asymptotic zone. Apart from these boundary conditions, q(p, z) is arbitrary. 

At a moment of time symmetry the constraint equations of General Relativity 
reduce to the single condition 

R = 0, (2) 

where R is the three-dimensional Ricci scalar. With Brill's choice of background 
metric, the Hamiltonian constraint (^) takes the form 

2 ( ip ( d 2 q t d 2 q 



which is solved for ip(p, z) once q(p, z) is given. In the current work we adopt the form 
of q(p, z) chosen by Miyama ||, 

q = ap 2 e- r2 (4) 

where r 2 = p 2 + z 2 , and the wave amplitude a is a free parameter. 

Since we require an axisymmetric, asymptotically flat solution which is reflection 
symmetric on the z = axis, we enforce the boundary conditions 

^(0,z) = 0, ^(p,0)=0, (5) 

on tp, together with a Robin condition at the outer boundary, 
dtp 1-ip 

~7T~ = , (6) 

or r 

which ensures that the hypersurface is asymptotically Euclidean. 

The solution for ip is calculated using these boundary conditions together with a 
second order accurate finite difference approximation to equation <M). 



N (number of grid points) 



(a) (b) 

Figure 1. (a) The conformal factor for a = 12 calculated using the Regge lattice 
with 600 X 600 vertices. The solution is shown on the first quadrant of the p — z 
plane, with the outer boundaries placed at p = z = 20. (b) The fractional 
difference, averaged pointwise across the domain, between the simplicial (ip r ) and 
continuum (ip c ) solutions for the conformal factor on an N X N grid. 



3. The Regge calculus model 

In this section we will briefly outline the axisymmetric tetrahedral lattice used to 
construct Regge initial data. Full details of the approach described here may be found 
elsewhere @. 

The simplicial equivalent of the initial value constraint at a moment of time 
symmetry, equation (0) , has been given by Wheeler , 

J2 L abtab = 0- (7) 

b 

This is the three-dimensional simplicial equivalent of R — at a vertex a, where L ao 
is the lattice edge joining vertex a to vertex b, and e ao is the deficit angle (product of 
the curvature on the edge by the area dual to the edge) about L ab . The summation 
is over all edges L ao which meet at vertex a. 

We mirror the continuum approach and perform a conformal decomposition on 
the lattice edges, since they are the simplicial equivalent of the metric. This gives 

L a b = Ipab ^ab 

where L a b is a freely chosen base edge length. It is most natural to define the simplicial 
conformal factor on the vertices of the lattice, such that ip a is the conformal factor 
at vertex a. The conformal factor ip ao acting on the edge L a & is then defined using a 
centred, second order approximation, 

iJab = ^ (V^a + i/fb) ■ (8) 

We assign the base edges L a b such that the lattice is aligned with the co-ordinate 
grids of the continuum approach. Figure || indicates how we subdivide each block 
in the co-ordinate grid to obtain a tetrahedral three geometry. Given the full set of 
physical edge lengths in this lattice, we calculate the deficit angles (curvatures) about 
each edge, and write down the simplicial initial value constraint, equation (fjj). This 
equation is then solved, together with suitable boundary conditions, for the conformal 
factor at each vertex of the lattice. 



(p,<|>,z) l ijk (p+Ap,(|),z) 



Figure 2. A portion of the simplicial three-lattice, which shows how each co- 
ordinate block is subdivided into six tetrahedra. After calculating the geometric 
quantities for this lattice, the limit A</> — > is taken, in order to obtain a 
discretization suitable for constructing axisymmetric initial data. 



To obtain a purely axisymmetric lattice, we take the limit as tends to zero, 
while demanding that each block is symmetric about (j> + \&4>- Retaining only the 
leading order terms of equation (f7j) in this limit yields the axisymmetric simplicial 
initial value equation at a moment of time symmetry. Full details of this calculation 
and the resulting equations may be found elsewhere Q . 

By aligning the lattice with a cylindrical polar co-ordinate system, we may use 
the form of the continuum base metric, equation (|l|), to specify the base edges in the 
lattice. Combining this with the conformal decomposition described above, and noting 
that after the limit has been taken there are only three non-zero edges per vertex, the 
physical lattice edges take the form 

*fc = [^e% +hk+k {^ + ^zl) 1/2 (9) 
hik= tyV].^ **k 

where all expressions are centred on the edge, and we have neglected the redundant 
index along the <j> axis. 

We are now able to calculate the simplicial initial data by solving equation (^) 
for the conformal factor ip. We use the same boundary conditions as before, applied 
along the p — and z — boundaries using power series expansions into the domain. 

The solution for i/j calculated using the Regge lattice, together with the boundary 
conditions described in section ^|, is shows in figure ^a. The solution displays excellent 
agreement with the previous calculations of Holz et al and Alcubierre et al [0, M, 
and the continuum model developed in section ^| Figure contains a rigorous 
convergence analysis which demonstrates that the Reggc solution converges to the 
numerical continuum solution as the second power of the lattice spacing. For further 
details, we refer the reader to Gentle Q|. 

Both solutions display good agreement with the previous calculations of both 
Holz et al {| and Alcubierre et al g . Table [l] contains estimates of the mass of the 
Brill waves, calculated using a least squares fit of the function 

M 

*a = l + i- (10) 

2r 
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to ip in the asymptotic zone. The mass M is found to converge rapidly with radius. 
The table also contains previous mass estimates |J, which again agree remarkably 
well with the current calculations. 

4. Trapped surfaces 

Our goal in the remainder of this paper is to locate any apparent horizon which exists 
in the Brill wave data. Apparent horizons are defined || as the outermost trapped 
surface on the spatial hypersurfacc. That is, the outermost closed, spacelike two- 
surface on which the divergence of outgoing null rays is zero. At a moment of time 
symmetry, this is equivalent to the outermost closed spacelike surface with extremal 
surface area. 

In order to construct the trapped surfaces we utilize the technique employed by 
Holz el al Q , where the search for an extremal two-surface was shown to be equivalent 
to the construction of geodesies in a certain two-dimensional space. 

Let us consider a path s(p, z) in the first quadrant of the p — z plane, which forms 
a closed surface when rotated about the z and p axes. The area of this surface of 
rotation is 

S = 27n/>Vp (dp 2 + dz 2 ) 1 ' 2 , (11) 

JPo 

and the criteria for a trapped surface is that this area is extremal, SS = 0. It is 
straightforward to show that this extremization is mathematically equivalent to finding 
the geodesies of the metric Q] 

ds 2 = 4tt We 29 {dp 2 + dz 2 ) , (12) 

subject to appropriate boundary conditions. If the geodesies are parameterized by A, 
the appropriate boundary conditions are 

p(X ) = and z(X ) = 0, (13) 
i(Ai) = and p(A x ) = 0, (14) 

where Ao and Ai are the endpoints of the curve, and p represents differentiation with 
respect to A. These boundary conditions ensure that the curves are continuously 
differentiable across the axes. 

The trapped surfaces are constructed using a shooting method. Given a value of 
p(\o), the geodesic is traced to Ai, where we check how closely the boundary condition 
is satisfied. Once a pair of p(\q) values have been found which bound the condition 
z(Xi) = 0, the bisection method is used to converge on the trapped surface. 

It would be interesting to construct an algorithm for finding trapped surfaces 
entirely within the framework of Regge calculus. An obvious approach at a moment 
of time symmetry would be to implement a direct search algorithm which minimizes 
the simplicial surface area of a topological two-sphere. For the purposes of the 
current work, we have assumed that the simplicial lattice corresponds to a section 
of a differentiable manifold, upon which we can apply the ideas outlined above. We 
use the conformal factor ip r , calculated using Regge calculus, to construct the geodesies 
of metric (|l2|), and hence construct trapped surfaces within the lattice. 
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Figure 3. (a) Trapped surfaces for a Brill wave with amplitude a = 12 at 
a moment of time symmetry. The outermost trapped surface is the apparent 
horizon. The calculations were performed using both the Regge (solid lines) 
and continuum (dashed lines) solutions on a 600 X 600 grid, and are in excellent 
agreement, (b) The point p(Arj) where the inner and outer trapped surfaces meet 
the p-axis is plotted as a function of wave amplitude, close to the critical amplitude 
a*. The quadratic fit suggests that a* £S 11.819 in the Regge solution (triangles), 
and a* pa 11.816 in the continuum initial data (squares), (c) Convergence of the 
trapped surfaces calculated using the two sets of initial data. For a Brill wave 
amplitude of a = 12, we show the fractional difference between the Regge and 
continuum values of p(Ao) as the grid resolution is increased. Both the inner and 
outer trapped surfaces are shown. The two solutions converge as the second power 
of the grid spacing. 



5. Apparent horizons 

Using the technique outlined in section ^, we now search for apparent horizons in the 
Brill wave initial data constructed in sections ^ and [| 

We find that in both the Regge and continuum solutions an apparent horizon first 
appears at a critical amplitude a*, which lies in the range a» e [11.81, 11.82]. As the 
amplitude a is increased above this critical value, two trapped surfaces are present, 
and these gradually separate. These results are in complete agreement with the recent 
calculations of Alcubierre et al M. 

Figure shows a projection of the trapped surfaces found for a wave amplitude 
of a = 12, calculated using both the continuum and Regge solutions. There is a slight 
difference between the surfaces calculated using these two independent sets of initial 
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data, but the discrepancy is found to reduce as the second power of the grid spacing. 
In figure ||b, we show the position p(\q) of both trapped surfaces along the p-axis as 
a function of the Brill wave amplitude. The quadratic fit to this data allows us to 
improve the estimate of the critical wave amplitude at which an apparent horizon first 
forms. We estimate that a* ~ 11.819 for the simplicial solution, and a* ~ 11.816 for 
the solution to the continuum equation. Both solutions predict that the single horizon 
at a = a* passes through p(Xo) ~ 2.291. 

The value of the critical amplitude a*, at which an apparent horizon first forms, 
differs from that obtained in an earlier calculation by some of us j|] . The discrepancy 
in the earlier work was due to the relatively weak convergence criteria used for the 
trapped surface finder. Qualitatively, the previous work agrees with both the present 
calculations and those of Alcubierre et al || . 

6. Future work 

We have used Regge calculus to confirm recent numerical findings regarding the 
formation of apparent horizons in a particular class of Brill wave spacetimes, while 
demonstrating that Regge calculus provides an alternative and competitive technique 
for use in numerical relativity. 

Brill wave spacetimes provide a challenging test-bed for the future development 
of simplicial gravity. Work is currently underway on the development of a simplicial 
trapped surface finder for Regge calculus, as well as the time evolution of the initial 
data constructed here, in both (2 + 1) and (3 + l)-dimensions. 
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